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Introduction 

In the upcoming era of permanent occupation of 
space by humans, the problem of adequate protec- 
tion from cosmic rays assumes greater and greater 
importance. An extensive research program is cur- 
rently underway (ref. 1) to address this concern. In 
particular, this program is concerned with the calcu- 
lation of cross sections describing the interaction of 
galactic cosmic ray nuclei with target nuclei. 

In the present work, we wish to show how to 
set up the general formulas for these cross sections 
(once the T-matrix elements are known) and how 
to obtain and handle the appropriate normalization 
volume elements which occur in the phase space and 
wave functions. More importantly, however, we wish 
to show how to calculate the correct phase space 
factors for two- and three-body final states. We shall 
present these calculations for both Lorentz- invariant 
and noninvariant phase spaces (ref. 2) and shall also 
show the equivalence of the two approaches that use 
energy derivatives and energy delta functions. 

The purpose of this report is to present the results 
of diverse phase space calculations in a consistent 
and unified manner. All too often the results of 
phase space calculations are presented in different 
ways depending on who is the author and what 
particular reaction is being studied. Thus, when one 
tries to calculate phase space factors for a particular 
reaction, it is often almost impossible to use the 
results of another author. Thus, a general formalism 
for the evaluations of phase space factors is needed 
that can be applied to any particular reaction. The 
specific results of other authors can be obtained 
through the use of our general results. 

Cross Section 

In scattering problems (as opposed to bound state 
problems) one introduces the concept of cross section 
< 7 , defined as (refs. 3 6) 



where w is the transition rate and is the incident 
flux defined as the number of particles per unit area 
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per unit time, that is, 

F t = riiPiV (2) 

where n, is the number of incident particles, p x is the 
incident probability density (i.e., per unit volume, so 
that n,p, is the number of incident particles per unit 
volume), and v is the incident velocity. Frauenfelder 
(ref. 3) points out that in most calculations the 
number of incident particles is normalized to one 
particle per volume V, that is, 

f ^*( r )^i(r) dV = 1 (3) 

Thus we shall take 


n i = 1 (4) 

Now in quantum mechanics, 


Pi = = y 

(5a) 

consistent with one particle per volume V , whereas 
in classical mechanics the probability is unity always 
so that again 

Pl y 

(5b) 

Thus, 


II 

(6) 

and the cross section can be written (ref. 2, 

p. 12) 

3 

^ I ^ 
II 
b 

(7) 

where the transition rate is given by 


w = — J |T fci | 2 dNips 

(8) 


where T ki is the appropriate matrix element and 
dNips represents the phase space elements. The units 
and volume normalization conventions are explained 
in appendixes A and B. 


The noninvariant phase space element (or density of states) is defined as (ref. 2) 


p N = cfNips {E,pi,P2, ■ ■ • i Pn) = 


V 


l(2nh)Z 


N—l / N 


N v / N \ / N 

,=1 7 V i=l 7 V »'= 1 


(9) 


(This definition of dNips differs slightly from that given on p. 388 of ref. 7). Frauenfelder (ref. 3) writes this as 


PN = 


JV -1 d 1 


L(2^)3J 




( 10 ) 


where E is the total energy of the final-state particles. In this version, the 6(E - has been replaced 

with d/dE and the 6 3 (p - E*iP.) has been taken into account by using in equation (10) rather 

than the n*i in equation (9). Although Frauenfelder’s version in equation (10) gives the same results as 
equation (9)*, we do not recommend its use because the inclusion of energy conservation is very obscure (being 
much more transparent and automatic in eq. (9)) and elementary errors are all too easy to make, such as 
bringing kinematic factors outside the momentum integrals when they should remain inside. We recommend 
using equation (9) to evaluate phase space factors, and we shall illustrate how to use it in the following sections. 

Two-Body Phase Space 

The two-body density is 

dNips(E; pi, p 2 ) = rf3 P2 6 ( E - £12 ) * 3 (p - E p * j 

= 1 535,/*, *(*-*») 

= ^^3 JJ Pi d Pi 6 ( E ~ E 12 ) 


with 

E\2 = E\ + E 2 


(12) 


From appendix C, we have 

in- (CM,,) (13) 

P\c* t, 1 + £/ 2 

which is valid in the center of mass of particles 1 and 2 (CM12). Note that E' 2 is identical to E 2 , but is now 
defined in a certain way as detailed in appendix C. Thus, equation (11) becomes 


dNips(£;pi,P2) 


V 


(27rfi) 3 




V 


(27T^)3 


hi 


£i4 JO 
p 'ENTe 2 da ' 


(14a) 

(14b) 


(Note the concerns raised in appendix D). However, there is something perplexing about equation (14b) We 
have supposedly integrated over d Pl ; yet Pl is still a variable in equations (14), both explicitly and also in 
E\ and E f 2 . This is resolved by considering the following. Energy conservation is required in equations (14) 
through S(E - E 12). (Momentum conservation is also included through <5 3 (p - Ei=iPi) ) In tbe frame, 
this implies that p\ = p\, as shown in appendix C. 
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Thus we need to express Pl as a function of E 12 , that is, 


Pi =Pi(E 12 ) 

(15a) 

and then because 


E i = E i (pi) 

(15b) 

E' 2 = E' 2 { P1 ) 

(15c) 

we have 


Ei = E X {E 12 ) 

(15d) 

e' 2 = E' 2 {E X2 ) 

(15e) 

l^l P fr^ZX he integrand ^ eqUati ° n (Ha) are fUnCti ° nS ° f * 12 ’ and the ~ 

^ 12 ) can then be 

Consider the nonrelativistic limit in the CM 12 frame. Then because p\ = p $, we have 


E u - p i + Pi Pi (cm t 

2 mi h 2m 2 ~ 2 P (CMl2 ) 

(16) 

so that 


Pl = Pl(£l2) = \f2 P E l2 

where P is the reduced mass of particles 1 and 2. Consequently, 

(17) 

P 1 — P\ __ PE U 
2m i m\ 

(18a) 

£> Pi _ PE X2 

2 m 2 m 2 

(18b) 

so that 


Pl E 1 ' + k-Jf^) 3/2 

(19) 

where M is the total mass of particles 1 and 2. Consequently, equation (14a) becomes 


«ps (E; P1 , K)m = ^^^ff {El2) W dEl2HB _ El2)dni 


(27 rh) 3 c 2 1 f dUl 

(20) 


An ^mnortanf 3 ^ “““! ^ "“'J “ <rgy and NR denoto -onrelativistic limit, 
dearly in equations (15)“ (Ml^Why ■* ,rf>lexin S (“•>) when things worked out so 
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This is illustrated by considering the relativistic evaluation of the two-body kinematic factor considered in 
appendix E. 

Three-body Phase Space 

The three-body density is 


dNips(E;pi,P2,P3) — 

/ 3 \ 

— / d 3 pi d 3 p2 d 3 P3 d{E ~ ^ 123 ) ( P ~ E Pl ) 

(2nh) b J v t=i ' 


- 

7-^rr c f d z pid 3 p 2 S{E - E 123 ) 
(27 m) D J 


= 

^ JJ pl dpi 6(E - E 123 ) dfl 1 d 3 p2 

(21) 

with 

E 123 = Ei + E 2 + E$ 

(22) 

From appendix F we have 



1 

I 

ElE 'z fCMioi) 

(23) 

Pl pic 2 ££ + (1 + P2/PI cos ^ 12)^1 


where 6 12 is the angle between particles 1 and 2. Thus equation (21) becomes 


V 2 

ffPl E l E Z *(p_ d 3 no 

(24) 

dNips(£;pi,P 2 ,P 3 ) - ( 271 -ft ) 6 

JJ c 2 E' 3 4- (1 + P 2 /P 1 cos0i 2 )E’i 

Evaluation of the kinematic factor gives 



Pl = Pl(-^123.P2) 

(25a) 

and then because 

Ei = E i(pi) 

(25b) 


E' 3 = E’ 3 (p i,p 2 ) 

(25c) 

we have 

El = Ei (Ei23> P 2 ) 

(25d) 


E 3 = E^(Ei23,P2) 

(25e) 

The p 2 that appears here is consistent with the integral over d 3 p 2 in equation (24). 
Evaluating the three-body phase space in the nonrelativistic limit, we have 



_ Pi pl . P1P2COS012 

123 “ 2pi3 4 2p 2 3 m 3 

(26) 

where 

m^rrij 

^ %3 TTlj -f- TTly 

(27) 
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Solving for p\ in equation (26), we have eliminated the variable p\ through 


Pi = Pl(-E’l23,P2) 


1 

P2 COS 012 

l P2 COS^l 2 \ 

M13 

m 3 \ 

V ™3 / 



~ 7 — W7. ^123 


(28) 


which is analogous to equations (17) and (E17). One now simply substitutes equation (28) into equations (25b) 
and (25c) and these results into equation (24). That is, one follows exactly the same procedure that we used 
for the two-body case in equations (18) through (20). Clearly the final result is an algebraic monstrosity and is 
best handled by computer code. One must be very careful about the signs for the square root in equation (28). 

We consider the relativistic three-body problem in appendix G. Ultrarelativist ic kinematics are presented 
in appendix H for both two-body and three-body problems. 


Dalitz Plots 

The Dalitz plot refers to three-body decays and tells us the kinematical distribution of particles 1 and 2 
due to phase space alone. Any deviation from the phase space diagram is due to variables in the T-matrix. An 
excellent intuitive discussion of Dalitz plots is given by Leon (ref. 4, pp. 145-147) and Pilkuhn (ref. 8, p. 160). 
The three-body density given in equation (21) is written now as 

dNips(E;pi,p 2 ,p 3 ) = ^— g J p\dpidn x J p| dp 2 d0 2 6(E - £ 123 ) (29) 

Recalling that #12 is the angle between particles 1 and 2, we can write (ref. 4, p. 57; ref. 9, p. 99; ref. 10, p. 329) 

V 2 r r 

dNips(£;pi,p 2 ,p 3 ) = jj pfp$d Pl dp 2 d(cos0i 2 ) 6{E - E 123 ) (30) 

In equation (21) we replaced dpi with a function of dE x 23 . In equation (30) we replace d(cos0 12 ) with a 
function of dE\ 22 > and replace dpi and dp 2 with functions of dEi and dE 2 , respectively. From appendix D, it 
is evident that 


However, equation (F4) gives 


E\ 

dp i = — dE\ 
Pi 

Eo 

dp 2 = — -dE% 
P 2 


G?(cos 012) = — ^ djEq23 

PlP2C l 


so that 


V 2 gyj-2 rr 

dNips(^ ; pi,p 2 ,p 3 ) = -J- Jj £i£ 2 £i23 dE x dE 2 dE 12 36 (E- E 123 ) 


(2ivh) Q c 2 

V 2 8 7T 2 


~ (2 j E ' E * dE 'dE2 


(31a) 

(31b) 


(32) 


(33) 


This, of course, is not the usual form of the Dalitz phase space, which is written using dLips where the product 
E\E 2 in equation (33) does not appear. (See the subsequent section on the Lor entz- invariant form of Dalitz 
plots.) 
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Recurrence Relations 

The density of iV-body final states can often be simplified in terms of products of two-body densities. 
Consider, for example, the reaction NN -+ N A -> NN 7r, where N represents a nucleon, tt a pion, and A 
an isobar. Normally we would simply write the three-body density dNips(E,p/v,PiV,P7r), especially if the 
square of the T-matrix | £ T kn G ni T ni \ 2 cannot be simplified. However, if the full T-matrix simplifies to 

n 

|T fcn | 2 |G n ^| 2 T nt | 2 , then the full cross section for NNtt can be written as a product of cross sections for 
forming NA and A — * Nir. Schematically this is 


°NN->NN tt ^ 


, TknG ni^ni 


(34) 


However, if 

a N N—*NN k * \Tkn\ 2 \Gni\ 2 \ T ni \ 2 ( 35 ) 

n 


then 


ff NN-*NNir ~ {&NN^NA){ l7 A—>Nn) 


(36) 


This equation shows intuitively that one might break the three-body factor dNips(E; p\ , p,\ : , Ptt) into a product 
of dWps(E;p N ,p A ) and dN\ps{E;p N ,pn). The result for three bodies is (ref. 2) 


dNips(£';pi,p2,P3) = d^ips(E;p c ,p^) dNips(E c \pi,p 2 ) dE c 


(37) 


and for four bodies is 

dNips(£;pi,p 2 ,P3,P4) = dNips(£; p c ,Pd) dNips(£ c ;pi,p2) dNips(E d ;p 3 ,p 4 ) dE c dE d (38) 

where dE c and dE d represent integration over the “nonobserved” particles. These equations can easily be 
generalized for N bodies. 

Finally we note that even though the process NN — *• NA — ► NNtt, where we have a well-defined 
intermediate state, may be occurring, the T-matrix may not break up as in equation (35), and so the best 
method would be not to use a recurrence relation, but rather the full three-body phase spaces developed earlier 
in the text. 

Lorentz-Invariant Phase Space 

We shall now repeat all the results thus far obtained and write them in terms of Lorentz- invariant phase 
space elements. (In this section we set h = c = 1.) The Lorentz- invariant phase space elements are defined as 
(refs. 2, 8, and 11) 

dUps(E-pi,p2, ■ • ■ ,Pn) = ( 27 r) 4 ~ 3JV '^ L ^ 4 (p ~ ( 39 ) 

This definition is (27r) 4 times the definition of reference 12. Note that the essential difference between dNips 
(eq. (9)) and dLips is the appearance of d' >> p l in dNips as opposed to d 3 p t / E l in dLips. This, in fact, makes the 
expressions for dLips somewhat simpler. 

Two-Body Phase Space 

The two-body density is 

dLips(E;p 1 ,p 2 ) = f d ^^6(E-E l2 )6\ P -p 12 ) (40) 
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with 


Pl2 = Pi + P2 

Equation (40) becomes 

rfLips(E;pi,p 2 ) = / ^^(^-^i 2 )dE 12 dni 

analogous to equation (14a). In the nonrelativistic limit this is 


dLips(£; pi , p 2 ) NR = /^i 


(41) 

(42) 

(43) 


analogous to equation (20). The relativistic expression is obtained by substituting equation (E17) into 
equation (42) and is much simpler than the rfNips expression. We give it as 

dLips(E;p 1 ,p 2 ) Rel = ^2 ^2 \JW 2 - (mi + m 2 )2] [£? _ (m , _ m2 )2] j ( 44 ) 

Three-Body Phase Space 

The three- body density is 

dLips(E; »„«,»)= ^ (45) 

(46) 


with 


P123 = Pi +P2 + P3 


Equation (45) becomes 


<*Lips(£; Pi , P2, p 3 ) 8(2t)5 J cosdl2)£l HE -E 123 ) dE 123 dfii d 3 p 2 (47) 

nonrZ!"? 11 ^ ‘ factor P roceeds b y substituting equation (28) into equation (47) for the 

nonrelativistic case. The relativistic case is considered in appendix G. 

Dalitz Plots 

The Dalitz plot form of the three-body density is treated in equations (29) through (33) for dNips The 
Lorentz- invariant form is evaluated in exactly the same manner. The result is particularly simple 

dLips(£; Pl ,p 2 ,p 3 ) = f dEi dE 2 ( 48 ) 

which is the usual form (refs. 4, 8, and 12) for the Dalitz plot. 

Recurrence Relations 

bodTes^ 0 ^ 2 mVafiant reCUrreDCe reIations are considered in references 2, 8, and 12. The result for three 


dLips(E;p 1 ,p 2 ,p 3 ) = — . dLips(E;p c ,p 3 )dUps(E c ;pi,p 2 )dEc 
and for four bodies (ref. 2) is 

dLips(£ , ,p 1 ,p 2 ,p 3 ,p 4 ) = ^ 27r )2 db ips(£; p c , p d ) dLips(E c ;pi, P 2 ) dLips(E d ;p 3 , p 4 ) dE c dE d 


(49) 

(50) 
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Concluding Remarks 

Methods of constructing phase space factors for two- and three-body densities have been presented for use 
in nuclear collision studies. Included in the discussion are methods for obtaining Lorentz-mvanant ones. The 
factoring of higher order phase space factors into products of lower order ones was also presented for use m 
simplifying nuclear cross section calculations. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
September 1, 1988 
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Appendix A 
Units 


Consider the units of the quantities in equa- where 
tions (5) through (8). The units of 7j., are seen from 


Tki = <fc|*|i> 

= f%t%dV (Al) 

and using box normalization as in 



(A2) 


we see that the volume terms in equation (Al) cancel 
so that the units of are Energy. (Note that 
instead of T*. t , we can have terms like 7fc n G m T n , 
where G is Green’s function with units (Energy) 
which cancel with the units of T m , so that T^ n G m T m 
has the same units as namely Energy). Looking 
at equation (9), we see that the units of dNips are 
always (Energy) -1 , no matter how many bodies are 
in the final state. Thus equation (8) can be written 
schematically as 


w oc 


^(Energy ) 2 g"— 


(A3) 


which gives the correct units for w of (Time) -1 . 
Clearly then <r, in equation (7), has units of 
(Length) 2 . 

In the above discussion, we have considered $ as 
a one- body wave function with normalization 1 j\/V . 
However, in a reaction 'L, is typically a two-body 
state (beam plus target) and an ,/V-body state. 
Consider a two-body final state. In second quantized 
form the (composite) two-body potential operator is 
(ref. 13) 


T = 2 EE <i ^i w>a « +fl i a l a k 

ij hi 

<^3 

k^l 


(A4) 


<ij\t\kl> = 



**( r l)*y(r2)«( r l.r2)^fc(r 1 )3- 1 (r 2 )c< 3 r 1 d 3 r 2 


(A5) 


so that obviously the units of <ij\t\kl> are still 
Energy. The operators a and a + have no units, 
so that T has units of Energy. When we take 
composite matrix elements as in nucleus- nucleus 
collisions, for example, <P 0 T 0 \T\P 0 T 0 >, we 
eventually obtain <P 0 T 0 \T\P 0 T 0 > in terms of 
< k p{P)h{T)\t\k p {P)l T {T)> (see eq. (36) of ref. 13) 
so that the units of <P 0 T 0 \T\P 0 T 0 > are still obvi- 
ously Energy. What happens when more than two 
particles are present in the final state as in the re- 
action NN -> N A -► NNn (ref. 14)? The matrix 
element is of the same form T kn G ni T ni , which can 
be written schematically as 


[<JV|(<7V|<7r|rtA>)|JV>] 

x [<JV|<A|G|7V>|7V>][<7V|< A|T|JV>|A r >] (A6) 

From our consideration of <ij\t\kl>, which we found 
to have Energy units, it is clear that G ni and T nl still 
have units of (Energy) -1 and Energy, respectively. 
The only concern is the extra 1 j\fV units carried by 
the newly created state <7r|. These units, however, 
are canceled by extra units carried in the pion- 
nucleon interaction. 


Let us summarize the ideas presented above. No 
matter how many bodies are present, T kl always has 
units of Energy. If T kl is replaced by terms such 
as ^ still has units of Energy, again no 

matter for how many bodies. The units of dNips 
are always Energy. Consequently w , as defined in 
equation (8), always has units of (Time) -1 and c 
units of (Length) 2 . 
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Appendix B 

Volume Normalization 


This appendix represents a continuation of the 
concerns raised in appendix A, where we established 
that our methods resulted in a consistent system 
of units. Here we are interested in the explicit 
cancellation of the volume factors V appearing in the 
expressions for o (eq. (7)), dNips (eq. (9)), and the 
wave function ^ (eq. (Al)). One might question why 
one should worry about volume factors V canceling 
when the units have already been established as 
correct. The answer lies in how one actually takes 
limits in box normalization. 

Consider the cross section expression from equa- 
tions (7) and (8) 


V2n f 

v h j 

\T ki \ 2 dNips 

(Bl) 

where (ref. 3) 



dNipsoc V N ~ l 

(For N > 2) 

(B2a) 

dNips cx V 

(For N = 1) 

(B2b) 

We write a schematically 

as 


o = {V K> 2 )(<fc| 2 V iV - 1 ) 

(B3) 

Now consider a one-body initial state and one-body 
final state. Then 

K > 2 

oc 1/V 

(B4a) 


<fc| 2 a 1/V (B4b) 


so that combining equations (B2b), (B4), and (B3), 
we see that the volume terms cancel. Now for a two- 
body state we do not write |i> 2 cx 1/V 2 because we 
use a two-body center- of- mass wave function with 
normalization \/\/V , so that equations (B4a) and 
(B4b) are still valid for the two-body center- of- mass 
reduced wave function. From equation (B2a), it is 
clear that the volume terms again cancel. 

Let us be a little more specific for the two-body 
case. Consider nucleus-nucleus scattering and write 

| %> = |$>iPa>|To> (B5) 

where represents the total initial state, rep- 
resents the wave function for the motion of the overall 
CM and is proportional to 1/ y/V ; and | Po> and |T 0 > 


represent the internal states of the projectile and tar- 
get, respectively, each with their own normalizations, 
say 1 /y/Vp and l/vA'T* Assume that the T-matrix 
separates as 

t = (B6) 

so that 


Now for example 

<p k \h\Po> « ^ p ~^f ktp ^° dVp (B8) 

so that all the V p volume terms cancel when the in- 
tegrals are performed (so that taking the limit as Vp 
approaches oo provides no difficulties). Let us con- 
clude then that all matrix elements connected with 
the internal states do not display V factors explic- 
itly. However, for the CM states in 
the cross section a is defined in terms of the normal- 
ization of which is why the V factor appears 

explicitly in equations (5) to (7). Clearly the l/VV 
factor of |$j> (when squared) always cancels the V 
factor of equation (7) because of the way they were 
defined originally. Similarly (in the two-body case), 
the l/VV factor of |$ fc > (when squared) cancels the 
V factor from the two-body dNips because again the 
dNips volume factor refers to the normalization of 

Now consider a two-body initial state with an 
iV-body final state. Clearly |$ t > volume factors 
cancel with the factor in equation (7) as before. The 
jV-body final state is written 

|fc> = 1 2-body CM>|7V-2 created particles> (B9) 

so that the final state has 2 + (N - 2), that is, N 
particles, in the final state. However, |2-body CM> 
carries only a 1/ \/V factor (as in eq. (B4b)), whereas 
the created particles each carry a 1 j\fV factor. (The 
two-body CM is like a single reduced particle as 
explained previously.) Thus, 

i / 1 \ N ~ 2 

(B10) 
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so that 

< fc ! 2 «^r (Bll) 

which always cancels the V N ~ X from dNips, as it 
should because of the way dNips is constructed. 

Finally we consider T kl replaced by terms like 
TknG n iT n i . The <k\ and |i> normalizations can- 
cel dNips and F u respectively, exactly as described 
above. Any extra 1/y/V terms arising should be can- 
celed by the integrals arising in the matrix elements 


exactly as described for equation (B7). We summa- 
rize as follows: the normalization of |t> (eq. (B4a)) 
must be written explicitly as 1/ \/V in order to cancel 
the normalization of Fp Similarly the normalization 
’ n ( e q- (Bll)) must be written explicitly to can- 
cel F 1 in dNips. All other normalizations (such 
as internal states or terms arising in T kn G ni T ni ) are 
canceled as in the integrals (eq. (B8)). The appar- 
ent inconsistency (i.e., |t> and <Jfc| not canceling in 
integrals) is due to our original definitions of F x and 
dNips in terms of | i> and |fc>. 


Appendix C 

Derivation of Equations (14) 

The total energy of particles 1 and 2 is 

E\2 = E\ + E 2 

— \/ (Pic) 2 + (m^ 2 ) 2 

+ y/ip^c^Tijn^c^ (Cl) 

However p } and p 2 are not independent variables 
because momentum is conserved. In the center-of- 
mass frame of particles 1 and 2, 



0 = Pl +P2 

(CM 12 ) 

(C2) 

so that 





Pi =P2 

(CM 12 ) 

(C3) 

Thus, 





^El2 _ dE i2 9 Ei + E' 0 
dpi dpi ~ Pl ° E x E' 2 (CMj 2 ) (C6) 

It is absolutely essential here to realize that we are 
now using 


E '2 = \/(Pi ~c)^~+Jm2C^ (C7) 

and not the form of E 2 in equation (Cl). That 
is, E 2 in equation (C6) is only a shorthand for 
equation (C7). We have emphasized this with the 
prime on E' 2 , which is identical to E 2 except that it 
is to be written as in equation (C7). 

Note also that equation (C6) is obtained from the 
Jacobians in appendix I (eqs. (15) and (16)) as 


so that E x 2 is a function of only p\ . that is, 
E12 = Ei 2 (pi) 

Thus, 


: 2 ) 2 

(C4) 

or 

dpi = 

dpi 

(C8) 

(C5) 


dpi = 

(dEi2\~ l 

(- 5 ?) 

(C9) 
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Appendix D 

Some Technical Points 


We raise some technical mathematical details now 
and illustrate them with reference to the two-body 
problem. They can be immediately generalized to 
the full N-body integral in equation (9). 

In equation (11) we had 

f d 3 p i d 3 P2 6 (E - Eu) S 3 (p - £ P* j 

= J d 3 pi 6(E — E\2) 

where the momentum-conserving delta function elim- 
inated the integral over d 3 p2 • Note that momentum 
conservation is thereby imposed on the right-hand 
side (RHS) of equation (Dl). 

Energy and momentum are related via 



Ex = \j{p\c) 2 + (raic 2 ) 2 

(D2) 

so that 

pi dpi = E\ dE\ 

(D3) 


Consequently 


J d 3 pi6(E - E 12 ) = JJp l dpi6(E - E l2 ) dn l 

PiE^iE -Ei 2 )dEi<K\i (D4) 


equation (D5) is radically different from the correct 
result in equations (14). The answer follows. 

In the left-hand side (LHS) of equation (Dl), all 
the variables are totally independent. In the RHS 
of equation (Dl), they are not independent, but 
are constrained by momentum conservation. If the 
energy delta function 

S{E-E 12 ) = 6{E-E l -E 2 ) (D6) 

were such that Ei and E 2 are completely indepen- 
dent variables, then the procedure in equation ( ) 

would be correct. However, E\ and E 2 are not inde- 
pendent; they are constrained by momentum conser- 
vation in equation (Dl). In fact in the center-of-mass 
frame 

6{E - Ei- E 2 ) =6 (e - ^/(pic) 2 -t- (mic 2 ) 2 

-\/(pi c ) 2 + (- 2 c 2 ) 2 ) ( D7 ) 


so that the combination 

J s(e~ \/(pic ) 2 + (mic 2 ) 2 

- ^/(pic) 2 + (m 2 c 2 ) 2 ^ dE\ ± 1 (D8) 



The subject of the present appendix is the following 
question: Why is it not correct to let 6{E - E n) 
eliminate the dE\ integral? Then we would obtain 

p\E\ S(E - £ 12 ) dE\ dQi « J piE\ dU\ (D5) 

After all, aren’t we just repeating the procedure used 
in equation (Dl)? 

This question is very subtle, but absolutely vital 
in obtaining correct results. The incorrect result in 


Thus the correct procedure (with reference to eq. (9)) 
is always 

1. Let the 6 3 (p - Dill Pi) freel y kiU a d * Pi 
integral. 

2. Then make sure that 6(E - ]Ci=l &*) an 
integral of the form d(£ili E t ) and not dE x . 

Forcing item 2 to occur always guarantees the 
correct kinematic factors. 
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Appendix E 

Two-Body Relativistic Kinematic Factor 

We define the metric tensor 


8 ^ = gT = 


10 0 On 


0-10 


0 


which transforms contravariant 4-vectors to 


0 0-10 

L0 0 0 -1 


covariant vectors A M via 


Ay. = gyi /A" 


Now 


A = 


rA°-\ 

A 1 

A 2 

La 3 J 


= (A 0 , A) 


so that 


Ay = 


0 

0 

0 

l — 


-A 0 ' 


’ A 0 - 

o-ioo 


A 1 


-A 1 

0 0-10 


A 2 


-A 2 

.0 0 0 -1. 


.A 3 . 


.—A 3 . 


and so (ref. 8, p. 2) 

Ay = (A 0 , -A) 

The scalar product of two 4-vectors is therefore (ref. 8, p. 2) 

A ■ B = A, ^ = A^Bv 
= A°B° - A B 

Another way to obtain the above result is simply to use complex Minkowski 

A = (A 0 , *A) 

but this has the disadvantage of hiding the true nature of the 


space 


The energy-momentum 4-vector 


is 


space-time metric (ref. 15). 


(El) 

(E2) 

(E3) 

(E4) 

(E5) 

(E6) 

(E7) 


Pn = (p°, -pc) = (E, -pc) 

or again we would use 


P = (E, ipc) 

The magnitude of all 4-vectors are Lorentz-invariant quantities, so that 


(E8a) 

(E8b) 


P = PP = p yP^ = E* - (|p| c ) 2 = (me 2 ) 2 


(E9) 


Now let us consider two-body kinematics. Aguilar-Benitez et al Cref 15 n isl a , , , 

body kinematics from the point of view of * Q ; no -l , , . ' ' 12, p ' 58 ) con sider two- and three- 

P °* V1CW ° f a Single P artlcle decaying into two or three bodies. They consider 


13 



the mattes in the test frame of the decaying partide ^hZXZ.2 

th^overaU (CM) ‘fame £TquiLen“ to the rest frame of the decaying part.de, we can take over 

th utt S -state partides f and , as 



Pl2 = ( e 12i -Pl2 c ) 

(ElOa) 

and the invariant mass is 

(M i2 c 2 ) 2 = p \ 2 

(ElOb) 

The overall CM frame is defined as 

pi2 = 0 = pi + P2 

(Ella) 

that is, 

P12 = (-£'12 . 0) — (A£i 2 c ,0) 

(El lb) 


Thus 


Pl2 = #12 = (^12C 2 ) 2 = (Pi + P2) 2 
-P 1 +P 2 + 2 Pl • P2 

= (mjc 2 ) 2 + (m 2 c 2 ) 2 + 2(E l E 2 - Pi ■ P2C 2 ) 


but in the CM frame 


Pl = -P2 


E 2 = E\2 — Ei 


so 


that 


Ei 2 = (mic 


2)2 + (m2C 2)2 + 2 [Ei(Ei 2 - El) + |pi| 2 ] 


= (mic 2 ) 2 + (m 2 c 2 ) 2 + 2 [^1(^12 - E l) + E 1 ( m l c2 ) 2 ] 


and thus 


E 2 x2 + (mic 2 ) 2 - (m 2 c 2 ) 
E l ~~ 2 Eu 


SO 


that (after some algebra) 


|Pll c = 


<J\El 2 - (m7+ m 2 ) 2 c 4 ] [E% 2 - (mi - m 2 ) 2 e 4 ] 


r2 

'12 


(E12) 


(E13) 

(E14) 


(E15) 


(E16) 


(E17) 


which is analogous to the result of reference 12 (p. 58) a 1th on eh the result is not as simple as the 

“ "“r^tJof momentnm or energy in the above anaiye. 
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Appendix F 

Derivation of Equation (23) 

The total energy of particles 1, 2, and 3 is 
-^123 = E\ -f E 2 + 

where, again, pi, p 2 , and p 3 are not independent variables because of conservation of momentum 
center-of-mass frame of particles 1, 2, and 3, 

0 = Pi + P2 + P3 (CM 123 ) 

so that 

P3 = Pi + v\ + 2piP2 cos 6 12 

where 0\2 is the relative angle between the momentum vectors of particles 1 and 2 
Thus 


^123 \J {P\ c )^ + (m^c 2 ) 2 -h \J (p 2 c) 2 4- (m 2 c 2 ) 2 

■h V (Pl c ) 2 + (P2 C ) 2 + 2piP2C 2 COS^i2 + (m 3 C 2 ) 2 
so that £ 123 is a function of both Pl and p 2 , but not p 3 , that is, 

^123 — Ei2s(pi,P2) 

Therefore, 

dE\23 _ d^!23 ^^123 dp 2 

dpi dpi dp A dpj 

but 

d P 2 

dpi 

because p 2 and p x are independent variables, and so 

dE \23 _ ggl23 Pl c 2 (Pi +P2cos0 12 )c 2 
rfpi <9pi ~ £1 + W 3 (CMi 23 ) 

where again it is essential to realize that 


E 3 = V (Pl c ) 2 + (P 2 c) 2 + 2p 1 p 2 c 2 cos0 12 + (m 3 c 2 ) 2 

Note again that 



using Jacobian notation of appendix I. 


(FI) 
In the total 

(F2) 

(F3) 

(F4) 

(F5) 

(F6) 

(F7) 

(F8) 

(F9) 

(F10) 
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Appendix G 

Three-Body Relativistic Kinematic Factor 


Our aim is to obtain the relativistic version of 
Pi ( £l23> P 2 ) gi ven i n equation (28). That is, we 
would like the three-body version of equations (E16) 
and (E17). The result is then simply inserted into 
equation (24). 

Proceeding as in appendix E, the overall CM 
frame is defined as 

p 123 = 0 = pi + P2 + P3 ( Gla ) 

that is, 

P123 = (■£'123)0) = (Mi23 c >0) (Gib) 

which are analogous to equations (Ella) and (El lb). 
Thus, 


so that 


,2x2 


^123 = ( m l c2 ) 2 + ( m 2 c2 ) 2 + ( m 3 c ) 

+ 2 (EjE 2 - |Pll|P2l c2 cos« 12 ) 
+ 2 [(Ei + E 2 ){Ei23 - Ei - E 2 ) 
+ (Pi +P2) 2c2 ] 

. 2)2 


(G4) 


= ( mjc 2 ) 2 + (m 2 c 2 ) 2 + (ms^) 

^E X E 2 - (m 1 c 2 ) 2 |p 2 |c 2 cos9i 2 


+ 2 1 


+ 2 


{El + E 2 ) [Ei 2 3 - E 1 - E 2 ) 


p 2 23 = E 2 23 = (Af 123 c 2 ) 2 = (PI + P2 + P3) 2 


+ E\ - (mjc 2 ) 2 + E\ - (m 2 c 2 ) 2 


= P? + P 2 + P 3 + 2pi ■ P2 + 2(P1 + P2) ' P3 

= (mjc 2 ) 2 + (m 2 c 2 ) 2 + (m 3 c 2 ) 2 + 2 (E, E 2 — pi P2^) 

+ 2 [(Ei + E 2 )E 3 - (pi + P 2 ) ' P3 c2 ] (G2) 

analogous to equation (E12), but in the overall CM 
frame 


P3 = -(Pl+P2) ( G3 ) 


2 \J~E% ~ ( m l c 2 ) 2 |p 2 1 cos $i 2 


(G5) 


analogous to equation (E15). 

Thus, by solving equation (G5) for E\ as a func- 
tion of £123 and |p 2 |, we have equations analogous 
to equations (E16) and (E17). Clearly the solutions 
are very complicated, and we shall not solve them 
explicitly here. 


Appendix H 

Ultrarelativistic Kinematics 


Let us evaluate the two- and three-body phase 
space factors in equations (14) and (24) using ultra- 
relativistic kinematics. Normally we have 


E 2 = (pc) 2 + (me 2 ) 2 

(HI) 

but in the ultrarelativistic limit 


mc 2 << pc 

(H2) 


so that 

E = pc (H3) 

is the ultrarelativistic expression relating energy and 
momentum (i.e., it is the expression for the photon). 

For the two-body final state, in the CM frame, 
the analogue to equation (16) is 

£l2 = (Pl c ) + (Pl c ) = 2 P lC 
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so that 


V 


Consequently 


Pl =Pl{E 12 ) = ~ 
Zc 


Ei = pi c = -E 12 


e 2 = -Ei 2 


so that 


(H5) 


(H6) 

(H7) 


rfNips(E; Pl ,p 2 ) UR = J Ef 2 6(E-E 12 )dE l2 dU 1 


V E 2 
(27rh)3 8c 3 


/ 


dfli 


(H9) 


For the three-body final state, in the overall CM 
frame, the analogue to equation (26) is 


E 123 = (p\c) + (p 2 c) + c\J p2 + p 2 + 2 pip 2 cos 0 12 

(H10) 

which we solve for pj as 


Pl 


El + E' 2 8c 


(H8) 


Equations (H4) to (H8) are the analogues to equa- 
tions (16) to (19). We substitute equation (H8) into 
equation (14a) so that the ultrarelativist ic expression 
for t/Nips is 


Pl — Pi(E\ 23 ,p 2 ) 

Ef 23 

2p 2 c 2 (£’i 23 + cos 0i 2 - 1) ^ Hl1 ^ 

analogous to equation (28) so that one would proceed 
in the same manner as described after 
equation (28). 


Appendix I 

Jacobian Transformations 


The change of variables from x and y to u and v 
where 


differentials via (ref. 9, p. 37) 


X = f(u,v) 


(11) 

dx\ dx 2 • • 

• dx n = 

d (xi,x 2 , 

d (yi,V 2 , 

•".*fi) j 

- vdy i 
* * ' j Vn) 

dy 2 • 

■ • . dy n 

II 

IT 


(12) 







(15) 

is accomplished through (ref. 16, p. 

556) 


with the Jacobian 






JJ <t>(z,y)dxdy = JJ 4>[f( 

u, t>), g(u, u)] dudv (13) 

o(u,v) K ' 








where the Jacobian is defined as 






dx\ 

<Jyi 

dx\ 

-&V2 

dx\ 

dyn 


d (x,y) 

dx dx 


(14) 

d(xi,x 2 ,- 

•■,X n ) 


dx 2 

(3yi 

pi ... 

oy 2 

dxo 


d(u, v) 

dy dy 

chi &v 


^(yi. 2 / 2 , • 

■' ,Vn) 





(A6) 

This may be generalized 

to transform any 

set of 





dxn 

^V2 ‘ * * 

pL 

oy n 
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Symbols 


A 

4-vector 

c 

speed of light, 3 x 10 8 m/sec 

d Nips 

noninvariant phase space element, 
MeV 1 

dLips 

Lorentz-invariant phase space element 

E 

energy, MeV 

^12i ^123 

defined in equations (12) and (22) 

fi 

incident flux, m~ 2 -sec“ 1 

gf tu 

metric tensor 


Green’s function, MeV -1 

h 

Planck’s constant, 6.58 x 10~ 22 MeV- 
sec 

|t > 

initial state vector 

|fc > 

final state vector 

TO 

mass, MeV/c 2 

M 

total mass, MeV/c 2 

n l 

number of incident particles 

N 

number of bodies in reaction; also 
nucleon 


P 

momentum vector, MeV/c 

r 

position vector, m 

t 

T-matrix operator 

T 

* ni 

T-matrix element 

V 

incident velocity, m-sec' 1 

V 

volume, m 3 

w 

transition rate, sec -1 

6 

Dirac delta function 

A 

isobar 

0 

angle, radians 

P 

reduced mass, MeV/c 2 

7 r 

pion 

Pi 

incident probability density, m - 

a 

cross section, millibarns 

4> 

one-body wave function, m~ 3 / 2 

$>, 9 

wave functions, m“ 3 / 2 

n 

solid angle, steradian 

* 

complex conjugate 
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